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1. Introduction 


The Electromagnetic Systems Branch (EV4) of the Avionic Systems Division at NASA Johnson 
Space Center in Houston, TX is studying the utility of surface acoustic wave (SAW) radio- 
frequency identification (RFID) tags for multiple wireless applications including detection, 
identification, tracking, and remote sensing of objects on the lunar surface, monitoring of 
environmental test facilities, structural shape and health monitoring, and nondestructive test and 
evaluation of assets. For all of these applications, it is anticipated that the system utilized to 
interrogate the SAW RFID tags may need to operate at fairly long range and in the presence of 
considerable multipath and multiple-access interference. Towards that end, EV4 is developing a 
prototype SAW RFID wireless interrogation system for use in such environments called the 
Passive Adaptive RFID Sensor Equipment (PARSEQ) system. The system utilizes a digitally 
beam-formed planar receiving antenna array to extend range and provide direction-of-arrival 
information coupled with an approximate maximum-likelihood signal processing algorithm to 
provide near-optimal estimation of both range and temperature. The system is capable of forming 
a large number of beams within the field of view and resolving the information from several tags 
within each beam. The combination of both spatial and waveform discrimination provides the 
capability to track and monitor telemetry from a large number of objects appearing 
simultaneously within the field of view of the receiving array. 

In this paper, we will consider the application of the PARSEQ system to the problem of 
simultaneous detection, identification, localization, and temperature estimation for multiple 
objects. We will summarize the overall design of the PARSEQ system and present a detailed 
description of the design and performance of the signal detection and estimation algorithms 
incorporated in the system. The system is currently configured only to measure temperature 
(jointly with range and tag ID), but future versions will be revised to measure parameters other 
than temperature as SAW tags capable of interfacing with external sensors become available. It 
is anticipated that the estimation of arbitrary parameters measured using SAW-based sensors will 
be based on techniques very similar to the joint range and temperature estimation techniques 
described in this paper. 

2. SAW RFID Background 

SAW RFID technology differs from more conventional integrated-circuit (IC) based RFID 
technology in many significant ways. For example, SAW RFID tags do not need to rectify 
incident electromagnetic power arriving at the tag. Instead, they modulate and re-radiate the 
interrogation signal directly using a series of reflectors printed on the surface of the device. The 
design and operation of a typical SAW RFID tag is illustrated in Figure 1 . 
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Figure 1. Design and operation of a typical SAW RFID tag. 

The principle benefit of this approach is that SAW RFID tags are capable of operating at much 
lower interrogator transmit power, and/or are capable of much greater range for a given transmit 
power than IC-based tags [1], In addition, SAW RFID tags have the capability to incorporate 
certain types of sensor telemetry along with the identification information. For example, the 
impulse response of a SAW RFID tag is highly sensitive to temperature. The temperature of the 
tag can be estimated directly by using a correlation operation to measure the time dilation (or 
contraction) of the impulse response. Finally, SAW RFID tags are very robust and can tolerate 
extreme conditions that would render IC-based tags inoperable. For example, they can withstand 
temperatures from cryogenic to over 310 C, they are inherently immune to ionizing radiation, 
and they have survived shock levels exceeding 1000 g. 

The temperature of a SAW RFID tag can be estimated by direct measurement of the time dilation 
(or contraction) of the tag impulse response. In particular, measurement of the time dilation of 
the impulse response at an arbitrary temperature relative to the response at a known reference 
temperature (usually 0° C) constitutes an observation of the temperature coefficient of delay 
(TCD) for the tag at its current temperature. Here, the term TCD refers to the mathematical 
function of temperature that quantifies the relationship between the relative time dilation of the 
tag response and the temperature of the tag, with respect to a fixed reference temperature. 
Although the TCD can theoretically be determined from the piezo-electric properties of the 
crystalline material used to manufacture the tag, it is more common (and probably more 
accurate) to estimate it experimentally. 

As an example, consider the waveform presented in Figure 2 below. This waveform represents 
the magnitude of a normalized version of the impulse response from a SAW RFID tag at a 
temperature of -73° C. The tag illustrated here is an example of the 40-bit Global SAW tag 
manufactured by RF SAW, Inc. [2]. There are a total of 22 pulses of slightly varying amplitude 
spaced intermittently over approximately a 1.6 ps time period after an initial delay of 
approximately 1 .2 ps. The exact spacing of the pulses, excluding the first and last pulse (which 
are used for synchronization) along with the phase of each pulse (which is not visible in Figure 
2) encodes the 40-bit ID of the individual tag. 
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Figure 2. Magnitude of normalized tag impulse response at -73° C. 

To illustrate the behavior of the RFID waveform a little more clearly, consider the waveforms 
illustrated in Figures 3 through 5. Figure 3 shows a detailed representation of one of the pulses 
from Figure 2 plotted showing both amplitude and phase information. As this figure indicates, at 
-73° C, each pulse lasts approximately 45 ns and oscillates at a frequency of approximately 2.44 
GHz. Similarly, Figure 4 shows a detailed representation of two pulses from the waveform 
depicted in Figure 2 shifted to baseband with a one-sided bandwidth of 80 MHz. In this case, the 
envelope of the signal is superimposed over the true signal waveform. This baseband equivalent 
waveform captures all of the information necessary to estimate time delay precisely from the 
pulses and eliminates the high-frequency oscillatory behavior that obscures some of the detail in 
the signal. Once again, the duration of each pulse is approximately 45 ns but the frequency of 
oscillation has been reduced to approximately 40 MHz. 
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Figure 3. Detail of pulse from Figure 2. 



Figure 4. Detail of two baseband pulses at -73° C. 
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Finally, Figure 5 shows the magnitude of all 22 pulses in the RFID response waveform at both - 
73° C (blue) and 21° C (green). It is clear from the figure that the higher temperature waveform 
has been dilated with respect to the lower temperature waveform. In fact, the dilation factor is 
approximately 1.0068, which corresponds to a change in dilation factor of approximately 

7.2x10 ^ per degree C. As it turns out, the TCD for these tags is very close to linear over large 

temperature ranges and 7.2x10 ~ represents a very good approximation to the slope of the 
corresponding line in the region from at least -100° C to 100° C. 



Figure 5. All 22 pulses of Global SAW tag impulse response magnitude at -73° C (blue) and 21° 

C (green). 

3. PARSEO System Description 

The hardware architecture for the PARSEQ system is illustrated in Figure 6. The components of 
the system are as follows. 

• National Instruments (NI) Chassis and Peripherals - The NI chassis contains a CPU 
board with USB and 100-Mbs Ethernet ports and a GPIB network interface board. The 
CPU board is used to control the operation of the overall system, perform all signal 
processing and beamforming operations, and provide an operator interface via an 
attached keyboard, mouse, and monitor. The GPIB network interface board supports 
communication between the CPU board and an RF switch matrix. 

• Agilent Network Analyzer (ENA) - The ENA is used to measure the frequency response 
of each of the 72 transmit-receive channels supported by the system. The ENA is 
operated in stepped-CW mode and measures the S12 channel response; that is, to 
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measure the frequency response at any particular frequency, a sinusoidal signal is 
transmitted continuously from Channel 1 of the ENA over a transmit antenna array while 
the amplitude and phase of the received sinusoidal response is measured on Channel 2 of 
the ENA through a single antenna element on a receive antenna array routed through the 
RF switch matrix. The ENA is connected to the NI Chassis via the Ethernet LAN, which 
carries both control and data traffic between the CPU board in the chassis and the ENA. 
The measurement of channel frequency response using the ENA is described in more 
detail in Appendix A of this paper. 

• Transmit Antenna Array - The transmit antenna array consists of 8 discrete circular patch 
antenna elements connected to the input via a 8-way power divider. The transmit antenna 
array is connected to channel 1 of the ENA via a low-loss, phase-coherent coaxial 
interconnect cable. 

• Receive Antenna Array — The receive antenna array consists of 72 discrete circular patch 
antenna elements. Each antenna element is individually connected to one of the input 
channels on the RF switch matrix via a low-loss, phase-coherent coaxial interconnect 
cable. 

• Keithley RF Switch Matrix - The RF switch matrix is used to connect individual elements 
of the receive antenna array to channel 2 of the ENA via low-loss, phase-coherent 
interconnect cables. The switch matrix is controlled by the CPU board in the NI chassis 
over the GPIB bus. 


r 



Tx Antenna 
Array 


National Instruments 
Chassis 

/\ 



Agilent Network 
Analyzer 



Qj (*) Q C-A 

O O G O 


QGQQGOQOO 


O O Q Q O 

o o o o 

J — — J 

OQQOQQQOO 

O Q O O & 

O G O O 


Q Q Q Q ^ Q Q Q (^) 


o o o o o 

G O O O 


OOOOGOOOO 


V 


Patch Antenna 



RF Switch 
Matrix 


Rx Antenna 
Array 


SAW RFID Tag 


Figure 6. Hardware architecture for EV4 PARSEQ SAW RFID interrogator system. 

The software architecture of the system, as currently configured, is illustrated in Figure 7. All 
software components of the system are hosted on the NI chassis CPU board. The software 
components of the system are as follows. 


7 








• LabViex v Script - The entire system is controlled by the LabView script. To acquire a 
single scan of data, the script sequentially sets the RF switch matrix to connect a single 
data stream from the receive array (i.e., a single antenna element) to Channel 2 of the 
ENA, commands the ENA to collect a sweep of the sampled frequency response for that 
channel, writes the frequency sweep data to a disk file, and then initiates the next step in 
the sequence. After collecting one full scan of data, the script passes control to an 
adaptive beamforming algorithm, which computes the desired beam coefficients from the 
sampled frequency-domain data and/or forms the desired beams using the computed 
beam coefficients. After forming all desired beams of frequency-domain data, the 
beainformer writes the beam data to a disk file and passes control back to the LabView 
script, which then passes control to a detection and estimation algorithm. The detection 
and estimation algorithm processes the beam data sequentially to estimate the range and 
temperature for each SAW RFID sensor tag detected in each individual beam. The range 
estimates along with the pre-determined azimuth and elevation associated with each 
beam are used to form a three-dimensional location estimate for each detected tag. After 
processing all beams, the detection and estimation algorithm writes the location and 
temperature estimates for all tags to a disk file and passes control back to the script, 
which displays the location/temperature estimates then continues the cycle by collecting a 
new scan of data. 

• Beamforming Algorithm — The beamformer operates in both an adaptive mode and a 
feed-forward mode. In the adaptive mode, the algorithm reads in the frequency-domain 
data for all of the receive array channels, computes and displays the spatial power 
distribution across the receive array, computes and stores the beam coefficients, azimuth, 
and elevation corresponding to a pre-selected number of the strongest beams, and 
displays the time-domain response for each selected beam. In feed-forward mode, the 
beamformer uses the previously computed and stored beam coefficients to form the 
desired beams in the frequency domain and writes the beam data to the disk file for use 
by the detection and estimation algorithm. In the current implementation, the beamformer 
also performs some preliminary signal conditioning operations prior to computing and 
storing the beam data. In future implementations, the signal conditioning operations will 
be moved to a separate signal conditioning block in order to facilitate performance tuning 
of the system. 

• Detection and Estimation Algorithm - The detection and estimation algorithm processes 
the beam-formed data and returns estimates of location and temperature for all tags 
detected in each of the beams. All processing occurs in the frequency domain, and tags 
are detected by correlating the received signal in each beam with a reference data base of 
tag templates measured at a pre-determined reference temperature. Only tags represented 
in the template data base can be detected by the current version of the software. It is 
hoped that future versions of the software will be modified to operate based on modeled 
rather than measured tag templates and to enable detection of all tags within the 
constellation of possible templates rather than a small subset of pre-determined templates. 


8 




Figure 7. Software architecture for EV4 PARSEQ SAW RFID interrogator 


4. Detection and Estimation 


The remainder of this paper is devoted to a fairly detailed description of the design and 
performance of the signal detection and estimation algorithm incorporated in the current 
PARSEQ system. An introduction to the basic concepts of detection and estimation for SAW 
RFID signals measured using a network analyzer is presented in Appendix A. Throughout the 
remainder of the paper, we adopt the notation used in Appendix A to describe the signal 
processing techniques. In particular, we assume that the received signal to be processed for each 
of the pre-selected beams takes the form R = S + N , where 

R = (Rq 9 R 1 ,...,R n _ 1 ), 

S - 

N = (N 0 ,N 1 ,...,N n _ 1 ). 


These three vectors represent the sampled frequency components for the total received signal in a 
particular beam, the received tag signal (which may be a composite of several tag response 
vectors), and the additive noise, respectively. The lowest sampled frequency, highest sampled 
frequency, and the number of sampled frequency components are all system configuration 
parameters. The noise samples are assumed to be complex-valued additive white Gaussian noise 

(AWGN) with mean zero and total variance a . As discussed in Appendix A, the variance is 
related to the power spectral density of the underlying real-valued continuous-time AWGN 
process observed at the network analyzer and the intermediate-frequency (IF) bandwidth setting 
on the analyzer. 
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The detection and estimation procedure described in Appendix A is a simple approach that works 
well in high signal-to-noise ratio (SNR) situations when there is little or no multiuser (i.e., 
multiple tag) or multipath interference within any beam. The approach implemented in the 
PARSEQ interrogator is considerably more sophisticated and performs much better in low SNR 
situations and in the presence of both multiuser and multipath interference. The implemented 
algorithm is a maximum-likelihood joint detection and estimation approach based on the wavelet 
transform. Although it is reasonably computationally intensive, it has been vectorized and 
optimized for fast computation using Matlab on a single processor. The algorithm can easily be 
partitioned for parallel implementation on multiple FPGA or DSP processors, and significant 
reductions in computation time can be realized in this manner. 

4.1. Basic Approach 

The algorithm operates by first computing a sampled version of the continuous-time wavelet 
transform of the so-called pre-envelope of the received signal with respect to the wavelets 
defined by the pre-envelopes of the RFID tag impulse response for each of the tags represented 
in the template data base. The wavelet transform [3] is precisely the right tool for our purposes 
because it provides a complete representation for an arbitrary time-domain signal in terms of 
time-shifted and time-scaled (i.e., compressed or dilated) versions of a particular template signal 
(the wavelet itself) 1 2 . Since we know that changes in temperature cause SAW impulse-response 
waveforms to scale in time, and changes in range cause the waveforms to shift in time at the 
receiver, the problem of interest is joint signal detection and parameter estimation for arbitrarily 
time-scaled and time-shifted versions of known template signals in the presence of AWGN. It 
follows from this that maximization of the wavelet transform of the received signal with respect 
to the wavelets defined by the SAW templates in the data base leads to a maximum-likelihood 
procedure for detection and range/temperature estimation of SAW RFID tag signals [4]. 

We assume that all tags are bandlimited to the interval [/ 0 ,/ 0 +R] Hz regardless of the 
temperature at which the tag is observed". If the tag frequency response is represented by the 
function //(/), then the impulse response takes the form 

*( 0 =£>(/)«“**. 

and the sampled tag frequency-response vector will take the form H = ( H 0 , T/j , . . . , H N _y ) , where 

H k =H{f 0 +j r k), k = 0,l,...,N-l. 

The corresponding pre-envelope of the tag impulse response is then given by 


1 This is analogous to classical time-frequency representations such as the short-time Fourier transfonn and the 
Wigner-Ville transfonns [2], but more mathematically rigorous. 

2 For our purposes, f 0 *2.36 GHz and B « 160 MHz. The interval chosen for processing is actually a system 
configuration parameter, but all tags currently available satisfy this constraint 
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h (t) = — X H k e 2 ^ /o ^' 

N k=0 


0 <t<T = N/B 


and (ignoring any time-domain aliasing resulting from sampling in the frequency domain) we 
have the relationship 

/7(f) = Re^/7 ( 7 )J . 

For the received signal vector R = (Rq, , . . . , R N _ 1 ) , we have the similar relationships: 

** =*(/„+■£*), k = 0,l,...,N-l, 


f (f ) = f- £ V ,2 ^ /o+S *^ , 0 <t<T = N/B, 

^ k = 0 

r(/) = Re[f(/)]. 

The wavelet transfonn of the pre-envelope of the received signal with respect to the wavelet 
defined by the pre-envelope of the tag impulse response (assuming it has zero DC component) is 
given by 

W r h (a, t) = — 1= f ?~(V)/7 dt, Va>0,Te'iH, 

Va J ~°° V a ; 

where the overbar indicates complex conjugation. If there are L tags in the template database 
with impulse-response signals {hi,h 2 ,...,h L } , and we know that exactly one tag is observed in a 

particular signal beam, then the maximum-likelihood estimate of the tag ID, temperature, and 
range is given by 

(/*,a*,x* j = argmax Re Wj ( a,x)l. 

|(/,cx,t)| 1=1,. . .,L\ a>0; xe9t j 

To perform practical joint detection and estimation using this approach without knowing a priori 
how many tags are present, we use a technique called successive interference cancellation (SIC) 
[5], which is described below. 

SIC, Stage 1. SIC is a two-stage joint detection and estimation process. To begin, we redefine 
the original pre-envelope received signal r(7) as the residual signal for Stage 1 processing and 

compute a sampled version of W r h , given by 

W r,/ ? (b J') = w r,h ( a / > T j ) > i = h 2, . . . , /; j = 1, 2, . . . , J , 
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for each of the impulse-response waveforms h e in the template database. The 

grid of sampling points z = 1,2,...,/; y = 1,2,..., jJ in the wavelet domain is arbitrary 

and can be chosen during system configuration. 

An example of this sampled wavelet transform is illustrated in Figure 8. For this example, the 
received signal is a time-scaled and time-shifted version of a single tag template embedded in 
AWGN at a signal-to-noise ratio (SNR) of 20 dB. In Part (a) of the figure, the chosen wavelet is 
also derived from the same template, while in Part (b), the wavelet was derived from a different 
tag template. In this case, the continuous wavelet transform was sampled at 800 time-shift points 
(approximately every 6.25 ns) and 41 time-scale points (approximately every 2.439° C). The 
time-scale and time-shift are referenced to the template itself; that is, the template is assumed to 
have been measured at a time shift of 0 and a time scale of 1 . As the figure indicates, when the 
chosen wavelet matches the received signal, the magnitude of the wavelet transform shows a 
strong peak at the sampled time-scale, time-shift point closest to the true time-shift and time- 

scale of the received signal, which in this case was the point ^0.998182, 0.667 128x1 0 -7 j. On 

the other hand, if the chosen wavelet does not match the received signal, the magnitude of the 
wavelet transform shows no strong peaks. 

Using the entire collection of sampled wavelet transforms for all tag templates in the data base, 
we can solve a discretized version of the desired optimization problem: 

(/*,/*,/*)= argmax Re 

{(/,/ j)|/=l,..., L\ 



using simple exhaustive search. If Re 


W rJ, f (<V ’ V ) 


does not exceed a pre-selected detection 

threshold, we declare that no tags are present in the current beam. If Re W r ; J+ J 

exceeds the detection threshold, we declare that some tag-like energy is present and perform 
secondary detection processing. Assuming tag /* satisfies the secondary detection criteria as 
well, then we declare that tag /* is present with the initial estimate of (a, t) given by j. 

To compute the final estimate of (a,x) for the detected tag, we solve the simplified continuous- 
time optimization problem 


(a*,x]= argmax Re 

»*. («*.*) 

{(a,x)|a>0; xe9?j 

L / J 


using a numerical optimization routine with starting value |ov*,x^j. The numerical 

optimization routine is the most computationally complex part of the algorithm, with a 
complexity of O ( TV 3 ] . 
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Figure 8. Sampled wavelet transform for received signal representing single tag response at SNR 

of 20 dB. 

SIC, Stage 2. If no tags are detected in the current residual signal during Stage 1, then detection 
and estimation processing for the current beam is terminated and processing proceeds to the next 
set of beam data. If a tag is detected in the current residual signal, then for the second stage of 
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processing, we assume that the estimated tag ID, temperature, and range from Stage 1 are all 
correct and reconstruct an estimated waveform corresponding to that particular tag signal by 
projecting the current residual signal onto the appropriately time-scaled and time-shifted version 
of the detected tag template from the data base. This reconstructed signal is then treated as a 
known interfering signal in the current residual signal, and a new residual signal is computed by 
subtracting the reconstructed signal from the current residual. This completes Stage 2 of the SIC 
processing, and control is passed back to Stage 1 to search for additional tags in the current 
beam. 

An example of the sampled wavelet transform before and after SIC processing is illustrated in 
Figure 9. For this example, the received signal is a superposition of time-scaled and time-shifted 
versions of two different tag templates embedded in AWGN. The SNR for each of the individual 
tag response signals, ignoring the interference due to the other tag, is once again set at 20 dB . 
The figure was computed using the wavelet defined by the tag template that produced the 
maximum observed sampled wavelet magnitude among the 1 1 templates in the data base. As 
expected, this was one of the two tag templates represented in the composite received signal. As 
before, the continuous wavelet transform was sampled at 800 time-shift points and 41 time-scale 
points, with the time-scale and time-shift referenced to the template defining the wavelet. Part (a) 
of the figure shows the sampled wavelet transform of the original composite signal during the 
first iteration of Stage 1 processing prior to any SIC processing. Part (b) of the figure shows the 
wavelet transform of the residual signal in the second iteration of Stage 1, after SIC processing 
during the first iteration of Stage 2 to remove the signal component corresponding to the first 
detected signal, which in this case was the signal matched to the chosen wavelet. As the figure 
indicates, prior to SIC processing the magnitude of the wavelet transform shows a strong peak at 
the sampled time-scale, time-shift point closest to the true time-shift and time-scale of the 

detected signal, which in this case was the point ^0. 998545, 0.800554 x 10 -7 J . After SIC 

processing, the magnitude of the wavelet transform shows no strong peaks, indicating that the 
signal component corresponding to the detected signal was successfully removed. 
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(a) Wavelet transform prior to SIC processing (first iteration of Stage 1). 



(b) Wavelet transform after SIC processing (second iteration of Stage 1). 

Figure 9. Sampled wavelet transform for received signal representing superposition of two tag 
responses at SNR of 20 dB using wavelet with strongest response (detected signal) during first 

iteration of Stage 1 . 
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For reference, we also include Figure 10, which again shows the sampled wavelet transform of 
the composite signal before and after SIC processing, but in this case, the figure was computed 
using the wavelet matched to the second of the two tag responses represented in the composite 
signal, which was the wavelet with the second strongest response in the first iteration of Stage 1. 
Part (a) shows the wavelet response during the first iteration of Stage 1 prior to any SIC 
processing to remove the first signal detected in Stage 1. The maximum magnitude of this 
response is slightly less than the maximum magnitude of the wavelet response illustrated in 
Figure 9(a), indicating that this tag response was not detected during the first iteration of Stage 1 
processing. Part (b) of Figure 10 shows the wavelet transform of the residual signal in the second 
iteration of Stage 1 , after SIC processing during the first iteration of Stage 2 to remove the first 
detected signal component. This response was in fact the strongest response among all the 
remaining templates in the data base (excluding the template for the first detected signal) during 
the second stage of Stage 1 processing. Hence, the second tag response represented in the 
composite received signal was detected during the second iteration of Stage 1. As the figure 
indicates, the magnitude of the wavelet transform still shows a strong peak at the sampled time- 
scale, time-shift point closest to the true time-shift and time-scale of the detected signal, which in 

this case was the point ^ 0.998 1 82, 0.667128 xlO -7 j . This indicates that very little of this signal 

component was removed during the first stage of SIC processing, which is the desired result. For 
this example, processing terminated during the third iteration of Stage 1, indicating that no 
further tag responses could be detected. 
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(a) Wavelet transform prior to SIC processing (first iteration of Stage 1). 
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(b) Wavelet transform after SIC processing (second iteration of Stage 1). 

Figure 10. Sampled wavelet transform for received signal representing supeiposition of two tag 
responses at SNR of 20 dB using wavelet with second strongest response during first iteration of 

Stage 1. 

4. Performance Evaluation 

In this section, we present an evaluation of the detection and estimation performance for the 
PARSEQ system. The performance of the system was evaluated both experimentally and using 
Monte Carlo computer simulations. The experimental performance was evaluated for the system 
as a whole, including the effects of the beamformer, but only estimation performance (as 
opposed to both detection and estimation performance) was quantified. The simulated 
performance was evaluated only for the detection and estimation algorithm itself, but both 
detection and estimation results were quantified. The relative efficiency of the estimation 
algorithm was evaluated with respect to the Cramer-Rao (C-R) lower bounds on the variance of 
unbiased estimates of both temperature and range. Analytical expressions for the (C-R) bounds 
were derived and evaluated numerically for comparison to the simulated performance results. 

4.1. Monte Carlo Performance Evaluation 

For the Monte Carlo performance evaluation of the detection and estimation algorithm, the tag 
template data base consisted of 1 1 tag waveforms that were measured experimentally by 
connecting each of the 11 reference tags to the SI port of the Agilent ENA network analyzer 
with an RF interconnect cable and recording the frequency response from each tag at 80 1 equally 
spaced points in the 160 MHz interval from 2.36 GHz to 2.52 GHz. The resulting measured 
frequency-response vectors were then processed to remove out-of-band noise, the two synch 
pulses, and other non-tag-related artifacts using the following procedure: 
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1. The tag response was filtered to the desired bandwidth by multiplying the measured 
frequency-response vector by a 801 -point Blackman frequency-domain window. 

2. The filtered frequency-response vector was converted into an over-sampled base-band 
impulse-response vector by padding the filtered frequency-response vector with zeros to 
a length of 65,536 samples and then performing an inverse FFT to produce a 65,536- 
point base-band impulse-response vector. 

3. The base-band impulse-response vector was windowed to remove the synch pulses and 
other non-tag-related transients by multiplying with a rectangular, unit-amplitude time- 
domain window with non-zero response in the interval from 1.26 ps to 2.63 ps and then 
normalized to have an average single- sample power of unity. 

4. The normalized impulse-response vector was converted back into the frequency domain 
by performing a 65,536-point FFT and then truncated back to 801 points to produce the 
final normalized frequency-response vector to be stored in the template database. 

The templates were all collected at room temperature, and the precise temperature of each tag 
during the measurement process was not recorded. Instead, an arbitrary reference temperature of 
20° C was assigned to each tag template for purposes of performance evaluation. Since the 
physical quantity being simulated and measured in this performance evaluation was actually time 
dilation of the tag waveform, from which temperature change was then inferred, the true 
temperature of each tag during the template measurement procedure was not relevant. An 
example of one of the resulting tag template impulse-response waveforms is illustrated in Figure 
11 . 
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To simulate the detection and estimation performance at a range of different temperatures and 
signal-to-noise ratios, each of the tags was “observed” at a different temperature and at many 
different SNRs by first numerically dilating or compressing the template impulse-response 

3 

waveform to the scale corresponding to the temperature assigned to that tag for the simulation 
and then adding simulated complex-valued, circularly symmetric, additive white Gaussian noise 
(AWGN) to the resulting tag frequency response vector. In particular, exactly one of the 
temperatures from the set {-5,0,5,10,15,20,25,30,35,40,45} was assigned to each of the 1 1 tags 

in the template database and all tags were simulated at 16 different SNR levels chosen at 1 dB 
increments in the range from 15 dB to 30 dB. For all simulations, the detection thresholds in the 
algorithm were fixed at levels selected experimentally to give good performance at a nominal 
SNR of 20 dB. Performance results for each SNR tested could have been improved by 
optimizing the detection thresholds for that SNR; however, such SNR-specific tuning is often not 
a realistic option for actual applications, so we used the more realistic approach of picking 
detection thresholds for a fixed nominal SNR value. 

4.1.1. Simulation Results 

For this set of simulations, a simulated received signal was derived by superimposing the 
simulated temperature-scaled frequency response vectors for multiple tags and then adding 
AWGN to produce the desired SNR level. That is, the SNR level was the same for all tags in the 
composite signal and was computed without consideration of the mutual interference produced 
by the multiplicity of tags. Detection and estimation results were simulated for composite 
received signals containing only a single tag, two tags, etc., up to the hill collection of 1 1 tags 
represented in the template database. The hill range of SNR levels from 15 to 30 dB was 
simulated for all tag multiplicity levels. 

The detection performance results for the Monte Carlo simulations are presented in Figures 12 
and 13 below. Figure 12 shows the probability of detection {Pd) for the full range of SNR values 
tested as a hinction of the number of total tag response signals present in the received signal. 
Here, Pd was computed by counting the number of tags in each simulated composite received 
signal that were correctly detected and dividing by the total number of tag response signals 
included in the composite. Notice that for the lowest SNR level of 15 dB, the detection 
probability was quite high (93.8%) when only a single tag was present but quite low (34.9%) 
when all 1 1 tags were present. However, as the SNR increases, detection probability grows 
monotonically regardless of the number of tags present, and for SNR levels greater than or equal 
to the nominal value of 20 dB, detection probability exceeds 94% uniformly across the entire 
range of tag multiplicity. In fact, at all SNR levels above 21 dB, detection probability uniformly 
exceeds 99%. 


3 ■ ■ • ^ 

For the simulation, tag TCD was assumed to vary linearly at the rate of 7.285714285714435 x 10 per degree C. 
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Figure 12. Probability of detection as a function of SNR and tag multiplicity. 

Figure 13 shows the probability of false alarm ( Pf a ) or false-alarm rate for the simulations. Here, 
Pf s was computed by counting the number of tags in each simulated composite received signal 
that were reported as detections but not actually present in the signal and dividing by the total 
number of tag templates not present in the composite. Notice that the false-alarm rate is less than 
5% unifonnly across the range of tag multiplicity for all SNR levels below 26 dB, and that it is 
less than 1% unifonnly across all SNR levels as long as the number of tags present in the signal 
is less than 5. However, for SNR levels above 25 dB, the false-alann rate begins to grow rapidly 
as the number of tags present in the received signal increases above 5, and for very high SNR 
levels, quickly reaches an unacceptable level. This is a result of the fact that some extraneous tag 
signal components always remain in the residual signal after the completion of the SIC procedure 
described above, and at sufficiently high SNR levels, the correlation between such extraneous 
components and the remaining tag templates can produce peaks in the wavelet transfonn that 
cannot be distinguished from true signal peaks. This problem can be ameliorated to some extent 
by introducing some adaptation into the selection of detection thresholds, but we are still 
experimenting with that procedure. 
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Figure 13. Probability of false alarm as a function of SNR and tag multiplicity. 

The estimation performance results for the Monte Carlo simulations are presented in Figures 14- 
17 below. Figures 14 and 15 show the mean squared error (MSE) for estimation of temperature 
and range, respectively, for the full range of SNR values tested as a function of the number of 
total tag response signals present in the received signal. As discussed in Section 4.1, both 
quantities were estimated simultaneously, but statistical analysis of the empirical results indicate 
that the two estimates are very weakly correlated in the low MSE regime. Also, while the 
evidence suggests that at least one of the estimators is not truly unbiased, which is discussed 
further below, the observed squared bias for both estimators was always at least an order of 
magnitude less than the observed MSE in the low MSE regime. Hence, the contribution of the 
bias to the empirical MSE of the estimates is insignificant in all scenarios where the temperature 
and range estimators perform well, so only the total MSE is presented here as a performance 
criterion. 

It should also be noted that estimation performance was tabulated only for correctly detected 
tags. That is, if a tag was not reported as detected, then no attempt was made to estimate either 
temperature or range, so it was automatically excluded from the statistical analysis. Similarly, if 
a tag was falsely reported as detected, and estimates of temperature and range were both 
computed, those estimates were excluded from the statistical analysis of estimation performance. 
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Figure 14. MSE for temperature estimation as a function of SNR and tag multiplicity. 
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Figure 15. MSE for range estimation as a function of SNR and tag multiplicity. 
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Notice that both figures indicate very good estimation performance across the entire range of 
SNR levels until the number of tags in the received signal exceeds 7. In fact, for all SNR levels 
in excess of 15 and tag multiplicities less than or equal to 7, the MSE for temperature is a small 
fraction of 1° C and the MSE for range is a small fraction of 1 meter. Notice also that for each 
SNR level below 26 dB, there is a very rapid increase in MSE for both temperature and range 
when the number of tags in the signal exceeds a certain level. Extrapolating just a bit, this would 
seem to indicate that there is a clear performance breakdown threshold for each SNR level, such 
that for tag multiplicities below that threshold, estimation performance is very good, and for tag 
multiplicities above that threshold, estimation performance quickly becomes very bad. For the 
lowest SNR tested (15 dB), this threshold seems to be 7 tags, while for SNR levels above 25 dB, 
the threshold seems to exceed 1 1 . 

We now consider the estimation performance in the absence of any tag mutual interference; that 
is, when only a single tag is present in the received signal and estimation performance is entirely 
a function of SNR. These results are presented in Figures 16 and 17, where the MSE for 
temperature and range are plotted, respectively, in comparison to the Cramer-Rao (C-R) lower 
bound on the variance (equivalently, MSE) of unbiased estimates for both quantities [6]. 
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Figure 16. MSE for temperature estimation as a function of SNR for a single tag response. 
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Figure 17. MSE for range estimation as a function of SNR for a single tag response. 

Two conclusions can be drawn from the results presented in Figures 16 and 17, namely, the 
estimates of temperature and range in the single-tag scenario are both very accurate and very 
efficient; however, some discussion is in order to clarify this. First, the actual MSE for both 
temperature and range is underestimated by the empirical MSE values in both figures at SNR 
levels 15 and 16. This is a result of the fact that not all tags in the simulated data set were 
correctly detected at these SNR levels. Although only a few tags were missed in each case, they 
were necessarily the weaker tags, so the average SNR of the detected tags was actually 
somewhat higher than the specified level. This implies that the empirical estimates of MSE in 
these two cases were necessarily biased below the true MSE of the estimators. If estimates of 
temperature and range for the missing tags had been computed and included in the statistical 
analysis, the observed values of MSE would have risen slightly. 

Second, the fact that the empirical MSE curve in Figure 16 actually falls below > the C-R lower 
bound curve is not an error. While, the C-R bound is a true lower bound on the MSE for any 
unbiased estimator of a parameter, it is possible in some cases for a biased estimator to have 
MSE lower than the C-R bound. Hence, the results in Figure 16 indicate that the bias of the 
temperature estimator studied here is not zero even though, as discussed above, it was observed 
to make an insignificant contribution to the error in the low MSE regime. Nevertheless, the MSE 
of the estimator is actually smaller than the MSE of an optimal (minimum MSE or MMSE) 
unbiased estimator, so the temperature estimator can be regarded as “ super-efficient ” relative to 
the MMSE unbiased temperature estimator. 

Finally, the C-R lower bound used for range estimation in Figure 17 is probably not the tightest 
possible lower bound for range estimation when temperature and range are estimated jointly, 
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since it actually represents the tightest possible C-R bound for estimation of range under the 
assumption that temperature is known a-priori, (i.e., temperature and range not estimated 
jointly). Even using this loose lower bound, the empirical MSE curve for the range estimates is 
uniformly less than 30% greater than the C-R lower bound curve across all SNR levels and 
approaches 20% at the highest SNR levels. Alternatively, shifting the empirical MSE curve to 
the left by 1 dB would put the curve uniformly below the C-R lower bound curve. Hence, while 
the range estimator does not outperform the MMSE unbiased range estimator, the SNR penalty is 
less than 1 dB relative to MMSE unbiased estimation of range alone. 

4.2. Experimental Performance Evaluation 

The experimental performance analyses performed on the PARSEQ system to date allow only a 
qualitative evaluation of the estimation performance of the system. They are included here in 
order to at least partially validate the simulation results as well as to convey some information 
about the estimation performance of the system in a dynamic environment. The discussion of the 
experimental performance results presented in this section is taken largely from an earlier 
technical report prepared by Mr. Reeve Ingle [7]. Minor editorial changes have been made here 
by the author. 

4.2.1. Experiment Description 

A single SAW RFID tag was placed in a temperature-controlled thermal chamber. An eight-foot 
length of Ml 7/1 58 coaxial cable was passed through an insulated port in the thermal chamber 
and used to connect the RFID tag inside the chamber to a single antenna element mounted on a 
tripod outside of the thermal chamber. This antenna was placed at a distance of two meters bore- 
site to the 72-element receiver antenna array. The transmit antenna was mounted on a tripod 
placed beside the receiver antenna array and pointed at the single-element antenna connected to 
the tag. The network analyzer IF bandwidth was set to 1 5 kHz during all thermal tests. 

The thermal chamber temperature was allowed to settle at 20° C to begin the experiment. A 15- 
minute thermal soak was performed at 20° C to allow the RFID tag to stabilize at this 
temperature. At this point, calibration parameters were set in the software to calibrate the system 
to generate a temperature reading of 20° C. Next, the thermal chamber temperature was 
decreased to -73° C at a constant temperature ramp rate (usually 2.5 °C/min). A thermal soak 
was again performed at -73° C for a minimum of 15 minutes to allow the tag to stabilize at that 
temperature. Subsequently, the thermal chamber temperature was increased back to 20° C at a 
constant temperature ramp rate (usually 4.5 °C/min). Throughout trials of this experiment, 
thermal chamber temperature data were time-stamped and recorded to an output file, and RFID 
tag range and temperature estimates were time-stamped and recorded to an output file. 

4.2.2. Results & Analysis 

As the thermal chamber temperature was decreased from 20° C to -73° C, the chamber was 
unable to maintain a constant temperature ramp rate over the entire temperature range. The 
temperature ramp rate decreased significantly as the chamber temperature approached -73° C, 
which is near the lower limit of temperature achievable with this chamber. Fortunately, the 
chamber was able to maintain a constant temperature ramp rate when the temperature was 
increased from -73° C to 20° C. Thus, estimation results are presented only for tests in which the 
temperature was increased from -73° C to 20° C. Also, due to a minor software bug that was 
later corrected, range measurements exhibited a temperature-dependent bias with a downward 
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trend as temperature increased. The actual range of the tag antenna was fixed at 2 meters 
throughout all trials. 

4.2.2. 1. Temperature Results 

Tests were performed at several different temperature ramp rates to test the ability of the 
PARSEQ system to track a time-varying temperature accurately. The results of these tests are 
illustrated in Figures 1 8-20 below. 



Time 


Figure 18. Comparison of measured temperature to actual chamber temperature for a temperature 

ramp rate of 4.5 °C/min. 

Figure 18 displays a plot of the actual thermal chamber temperature (red trace) and the 
temperature measurements from the RFID tag (blue trace) as the temperature was increased from 
-73° C to 20° C at a temperature ramp rate of 4.5 °C/min. As the figure indicates, the RFID 
temperature measurements began approximately 0.2° C above the actual chamber temperature 
and ended approximately 1° C above the actual chamber temperature after a 25 -minute thermal 
soak at 20° C. Note that the system tracked the actual temperature measured inside the chamber 
to within an accuracy of less than 1° C over the entire ramp but overshot slightly at the end of the 
ramp. 

It should be noted that switching between all 72 elements in the receiver antenna array requires 
approximately 20 seconds, and the thermal chamber temperature increases by approximately 1.5° 
C during that time at a ramp rate of 4.5 °C/min. Thus, each antenna element receives a tag 
waveform that is scaled differently than the tag waveform received by the previous element, and 
the output of the beam-forming program is “smeared out” in time. The output temperature 
computed by the detection and estimation algorithm based on the beam-formed data seems to 
equal roughly the average chamber temperature during the 72-element sweep. During subsequent 
trials, the temperature ramp rate was increased to try to determine the maximum rate that the 
system is capable of tracking. That is, above a certain temperature ramp rate, it is to be expected 
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that the beam-formed signal will be so badly distorted that the estimation algorithm will no 
longer be able to either detect the tag or produce an accurate temperature estimate. 
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Figure 19. Comparison of measured temperature to actual chamber temperature for a temperature 

ramp rate of 6 °C/min. 
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Figure 20. Comparison of measured temperature to actual chamber temperature for a temperature 

ramp rate of 1 1 °C/min. 

Figure 19 shows the thermal chamber temperatures and RFID tag temperatures after the 
temperature ramp rate was increased to 6 °C/min. At this rate, the thermal chamber temperature 
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changes 2° C during the time it takes to sweep all 72 antenna elements. The RFID temperature 
measurements generally overstated the actual temperature in the chamber slightly but remained 
within approximately 2° C of the actual thermal chamber temperature throughout the test. 
Similarly, Figure 20 shows the results of the experiment when the temperature inside the 
chamber was increased at the maximum possible rate, which was approximately 1 1 °C/min at the 
peak. At this rate, the thermal chamber temperature changes approximately 3.7° C during the 
time it takes to sweep all 72 antenna elements. Once again, the temperature measurements 
overstated the actual chamber temperature slightly but remained within approximately 2° C of 
the actual temperature throughout the test. 

Hence, although it is clear that the system will fail to track accurately at some rate of temperature 
change, the PARSEQ system was able to track temperature reasonably accurately at all ramp 
rates tested in these experiments. 

4.2.2.2. Range Results 

An example of one set of range measurement results is illustrated in Figure 21 below. These 
measurements were taken as the chamber temperature was increased from -73° C to 20° C at a 
ramp rate of 4.5 °C/min. As discussed above, a software bug that has since been corrected 
created a temperature-dependent bias in the range estimates that makes it difficult to draw many 
conclusions from these results. Nevertheless, if the linear trend is removed from the displayed 
data, it is clear that the variability of the range estimates falls within a range of approximately 
±0.2 meters over the entire temperature ramp. 
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Figure 21. Plot of range estimates vs. temperature for a temperature ramp rate of 4.5 °C/min and 

fixed actual range of 2 m. 

5. Conclusion 

This paper briefly summarizes the design of the NASA JSC PARSEQ long-range SAW RFID 
interrogation system and presents a detailed description of the design and performance of the 
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signal detection and estimation algorithms incorporated in the system. Numerical performance 
results indicate that the detection and estimation algorithm can detect and discriminate between 
at least 11 superimposed tag waveforms with very high probability in composite signals 
corrupted with AWGN at SNR levels as low as 20 dB. At the same SNR levels, the algorithm 
produces extremely accurate estimates of both temperature and range for up to 8 tags 
simultaneously, and continues to produce reasonably accurate estimates with as many as 1 1 tags 
present simultaneously. For all SNR levels above 25 dB, temperature and range estimates are 
extremely accurate for all tag multiplicity levels tested. Further, for signals containing only a 
single tag, the temperature estimation outperforms the best possible unbiased estimator (in terms 
of MSE) for the entire range of SNR levels tested, while the range estimation performance comes 
within 30% of that of the best possible unbiased estimator over the same SNR range. 

In terms of experimental performance, the available results qualitatively validate the numerical 
performance results and indicate that the PARSEQ system can reliably estimate temperature and 
range in real time and track temperature variation at rates up to 1 1 °C/min. 
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Appendix A 

Notes on Network Analyzer Detection and Estimation for SAW RFID Applications 


A.l. Introduction 


Assume that a system has real-valued impulse response /?(f) with corresponding frequency 
response //(/) that satisfy 

h(t) = 0, \/t < 0 and t > — , 

B 


and 

H{f) ~ 0 for all |/| < / 0 and \f\>f 0 + B, 

where B « / 0 is the system bandwidth and N is the desired number of frequency samples 
(bins). These conditions imply that h(t) can be recovered with little or no error from the 

frequency- domain samples //(/ 0 +y k) for k = 0, 1, 2, [N- 1) . Furthermore if we have the 
system input- output relationship 
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Also, using a network analyzer to sample the system frequency response as we are doing, the 
output from the system can be represented in the form 
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for k = 0, 1, 2, (#-1), where S k = + and the N k are i.i.d. CA(0,223jpAf 0 ) • 

A.2. Signal-to-Noise Ratio 

Assuming we are trying to detect the true unit impulse response /?(/) in the time domain in the 
presence of additive white Gaussian noise (AWGN) with power spectral density jV 0 /2, the 
signal energy is given by 

*=J „>(0f*=fZ4 2 . 


k= 0 


and the optimal detector performance for this signal is determined by the signal-to-noise ratio 
(SNR) 


2 B ^ .2 
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If we normalize this assuming that the system is excited with a band-limited impulse having 
energy £ rather than a perfect unit impulse function, which has infinite energy, then the signal 
energy is given by 

c N - 1 

c ^ .2 


and the normalized version of SNR, becomes 
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Similarly, assuming the analyzer output Fourier coefficients are known, the perfonnance of 
the optimal detector based on the frequency-domain sampled signal vector R = S + N , where 

R = 
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is determined by the alternative SNR 

IISII 2 1 N-l 
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If we again normalize this assuming that the system is excited with time-limited sinusoids having 
total energy £ rather than pure sinusoids having infinite energy, then the normalized version of 
SNR 2 becomes 


SNR 2 = 


N - 1 


N'N 0 N-N 0 ^ 0 


X4 


which is identical to SNRi 


Finally, the performance of a simple energy detector based on the signal vector R = S + N is 
loosely 4 related to the so-called deflection ratio , which is given by 


SNR 3 = 
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which can again be nonualized to give 
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A.3. Simple Estimation of individual RFID Tag Delays 

Suppose that we measure the frequency response of a SAW RFID tag and we wish to estimate 
the individual pulse delays in the tag response directly from the frequency response. In this case 
the frequency response of interest is assumed to take the form 


H k = h (/ 0 +jik) = y * =o, i, 2, .... (jv-i), 

j = o 


4 For a complete discussion of the problem of SNR and error probability for energy detectors, see “On Generalized 
Signal-to-Noise Ratios in Quadratic Detection,” R. J. Barton and H. V. Poor, Mathematics of Control, Signals, and 
Systems , 1992. Generally, if N is large and the Central Limit Theorem holds, the deflection ratio corresponds to the 
usual SNR for this problem and is closely related to probability of error. 
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where x ■ , j = 0, 1, . . . , P - 1 , represent the unknown delays of the individual reflectors on the tag, 

and Aje l0j , j = 0, 1, . . . , P - 1 , represent the known reflection coefficients of the reflectors. Note 

that with this approach, we are implicitly assuming that the reflection coefficients are 
independent of frequency across the entire bandwidth of the device. Changing this assumption to 
make the algorithm more accurate is straightforward but complicates the notation and is often 
neglected in practice. 

To estimate the sequence of delays in a very straightforward manner, we start be reconstructing a 
baseband, analytic version of the impulse response given by 
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It follows that the relationship between hit ) and is given by 

h(t) = Re h[t^e 2l ^ ot = v(7)cos(27t/ 0 7)-y(/ t )sin(27i/ 0 / t ). 


where x(7) = Re /?(/) andy(7) = Im /7(f) . 


Note that if the delays satisfy the condition 
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Furthermore, if f = x . + A where |a| < 1/5 , then 
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It follows that the delays x j = 0, 1, P-1 , can be estimated by finding all of the points 
tj=Tj+Aj where the observed (i.e., noisy) version of the function | h(t) has a relative 
maximum above some detection threshold and then solving for A • using the relationship 

Hence, the final estimate for each x ■ , j = 0, 1, . . . , P — 1 , is given by t ; = tj - A ■ , where 


A y = 


arg|^/7 j represents the phase of j , and the expression [jc] 2jt represents the number jc 
reduced modulo 2 ti. 

In order to evaluate h(t ) at enough values of t to give high-resolution estimates of the 
significant relative maxima, we use the inverse FFT. We start with the vector 
H = ( H 0 , Hi,..., H n _ i ) and pad with zeroes to a length of M > ( 2 f 0 N )/5 , where M is a power 

of 2. Let H = . . , H M _ X ) , where 
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Hence, the vector h = [h 0 , h x , ..., h M _ x j , where h m - j for m = 0, 1, . . M - 1 , is a scaled 

version of the inverse FFT of the vector H = ( H 0 , H x ,..., Hm -\ ) • 

A.4. Estimation of Temperature From Tag Delays 

One approach to estimating the temperature of a SAW RFID tag, although not a very good one, 
is to start by measuring the individual pulse delays in the tag response for some subset of tag 
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pulses. The temperature of a SAW RFID tag can be estimated by measuring the tag delays if the 
temperature coefficient of delay (TCD) for the tag is known relative to a fixed reference 
temperature and the tag delays are also known at the reference temperature. The TCD of a tag is 
a function that represents the relative change in delay for any fixed reflector on the tag, with 
respect to the delay at the reference temperature, as a function of the change in temperature from 
the reference value. It can be represented equivalently as the relative change in the acoustic wave 
propagation velocity along the tag as a function of change in temperature, but the 
characterization we give here seems to be the one used most commonly in the literature. 

If we let x(T) represent the delay associated with a fixed reflector at an arbitrary temperature T, 
and we let T = T 0 represent the reference temperature for the device, then the TCD can be 
approximated for most purposes over fairly large temperature ranges as 

ldlfiM*A(T-T 0 ) + B{T-T 0 ) 2 , 

x ( T o) 

where A and B are constants that must be determined experimentally for the particular class of 
devices in question at the reference temperature. Letting Tg=T(r 0 ), At = x(r)-x(r 0 ) , and 

AT = T —Tq, this becomes 

f(AT) = — ~A(AT) + B(AT) 2 . 

T o 

It follows that for any known set of reference delays x 7 - 0 = x 7 - (T 0 ) , i = 0, 1, ..., P — 1 , we must 
have 

Ax ? 

-t- = f(AT)*A(AT) + B(AT) , Vi = 0, 1, P-1 . 

T i,0 


Hence, to estimate the current temperature on a tag, one approach is to first use the algorithm 
described in the previous section to estimate the current value T 7 -=x 7 (r) for each delay, 

compute the corresponding estimated change in delay At 7 - = x 7 - — x 7 0 , and then solve the 
minimization problem 


min 

AT 



At, 
t i, o 


a(at)-b(at) 2 


to find the best least-squares estimate of AT . 

It is easy to see that there are tw o possible solutions to this problem, given by 



where 
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Ax = — 2 — 

p h t »,o 

is the average relative change in delay. Evaluating the mean squared error (MSE) for these two 
solutions gives 


p- 1 

2 “ 

2 p- 1 

_ 

MSE* = X 

7=0 


= Z 

7=0 

— -At 

T i, 0 


in both cases. Hopefully, only one of these two solutions makes sense physically. 
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